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We explore how the gravitational self force (or "radiation reaction" force), acting on a pointlike 
test particle in curved spacetime, is modified in a gauge transformation. We derive the general 
transformation law, describing the change in the self force in terms of the infinitesimal displacement 
vector associated with the gauge transformation. Based on this transformation law, we extend the 
regularization prescription by Mino et al. and Quinn and Wald (originally formulated within the 
harmonic gauge) to an arbitrary gauge. Then we extend the method of mode-sum regularization 
(which provides a practical means for calculating the regularized self force and was recently applied 
to the harmonic-gauge gravitational self force) to an arbitrary gauge. We find that the regularization 
parameters involved in this method are gauge-independent. We also explore the gauge transforma- 
tion of the self force from the harmonic gauge to the Regge- Wheeler gauge and to the radiation 
gauge, focusing attention on the regularity of these gauge transformations. We conclude that the 
transformation of the self force to the Regge- Wheeler gauge in Schwarzschild spacetime is regular for 
radial orbits and irregular otherwise, whereas the transformation to the radiation gauge is irregular 
^ ' for all orbits. 

hj ; 

\Q ■ 04.25.-g, 04.30.Db, 04.70.Bw 



> ■ I. INTRODUCTION 

\o 

( -—i ' Recent works, by Mino, Sasaki, and Tanaka Q, and by Quinn and Wald || (MSTQW), established a formal 

framework for calculating the local gravitational self force acting on a pointlike particle in curved spacetime. In these 
^^ ■ works, a particle of small mass m was considered, whose gravitational field may be treated as a small perturbation to 
the (vacuum) background metric. Such a finite-mass particle does not follow a geodesic of the background geometry, 
as its interaction with its own gravitational field gives rise to the exertion of a "self force" . In the above works, a 
general formal expression was obtained for the 0(m) self-force correction to the geodesic equation of motion. 

From the astrophysical point of view, the pointlike particle model and the self-force phenomenon may be applicable 
to binary systems with an extreme mass ratio. Of particular relevance are binary systems composed of a solar-mass 
(3J), compact object orbiting a supermassive black hole (of the kind now believed to reside in the cores of many galaxies). 
Such systems are expected to serve as main targets for the proposed space-based gravitational wave detector LISA 
(the Laser Interferometer Space Antenna), specializing in the low frequency range below lHzl. Knowing the local 
self force would be necessary, in general, for describing the orbital evolution in such systems, and, eventually, for 
?H \ characterizing the consequent waveform of the gravitational radiation emitted. 

When considering a model of a pointlike particle, one unavoidably encounters divergent quantities: the perturbed 
metric diverges at the location of the particle, and the "bare" self force associated with the metric perturbation turns 
out indefinite. One then has to deal with the fundamental issue of regularization; namely, extracting the correct, 
physical self force from the (indefinite) expression for the bare self force. The combined works by MSTQW present 
three different physically-motivated methods of regularization, all yielding the same formal expression for the physical 
self force F s " lf . This expression can be written in the schematic formFI 

pa _ pa pa fi\ 

self bare inst ' \ / 

Here F^ aic is the "bare" force, derived by applying a certain differential operator [see Eq. (]2j) below] to the full 
metric perturbation produced by the particle, and F? st is the singular piece to be removed. According to MSTQW 



1 Strictly speaking, both quantities on the right-hand side of Eq. (hi) are indefinite as they stand. In practice, one actually 
defines these two quantities as vector fields in the neighborhood of the particle. Then, the self force F£, l{ is obtained by taking 
the (well defined and finite) limit of the difference jr a ( baro ) _ i? Q ( 1Ilst ) as the particle is approached. For simplicity, we shall not 
use here this more strict formulation. 



analyses, this singular piece is to be constructed from the local, "instantaneous" part of the metric perturbation in the 
harmonic gauge, i.e., the part directly propagated along the light cone. The finite difference F" gxe — F^. represents the 
effect of the "tail" part of the particle's gravitational perturbation — the part scattered off spacetime curvature before 
interacting back with the particle. [The result by MSTQW is formulated in terms of the retarded Green's function. 
The bare force is then expressed as an integral (of a certain combination of Green's function derivatives) along the 
entire worldline of the particle, while the instantaneous part F^ st arises from integration along an infinitesimal, local 
piece of the worldline, that contains the momentary particle's location.] 

The first direct implementation of MSTQW's prescription for an actual calculation of the self force was carried out 
recently by Pfenning and Poission B , who considered the motion of a particle in a weakly curved region of spacetime 
(Pfenning and Poission also calculated the electromagnetic and scalar self forces acting on a particle endowed with 
electric or scalar charges, respectively). To allow calculation of the gravitational self force in strong field as well, 
Barack || recently introduced a method of multipole mode decomposition, based on the formal result by MSTQW. 
This method of "mode sum regularization" was previously developed || and tested J7J for the toy model of the scalar 
self force. We comment that a different mode-sum approach to the gravitational self force was proposed by Lousto 

The gravitational self force — unlike its electromagnetic or scalar counterparts — is a gauge-dependent entity. This 
statement means that the value of the self force is changed, in general, when the metric perturbation to which it 
corresponds is being subject to a gauge transformation (i.e., an infinitesimal coordinate transformation). If fact, the 
self force can be nullified along any segment of the worldline by a suitable choice of the gauge. Thus, any expression 
for the self force would be meaningless, unless one is provided with the information about the gauge to which this 
force corresponds. In MSTQW's analysis, the construction of the self force is formulated within the harmonic gauge, 
and the resulting expression (fil) therefore describes the harmonic gauge self force. Likewise, all implementations of 
MSTQW's analysis considered so far j|,|[ have been confined to the framework of the harmonic gauge, and have 
yielded the harmonic gauge self force. 

It is of great importance to understand the gauge dependence of the self force and to figure out how to construct 
it in gauges other than the harmonic: From the theoretical point of view, characterization of the self-force's gauge 
dependence is essential for a better understanding of the self force phenomenon; From the practical point of view, 
the harmonic gauge is not the most convenient one for actual calculations, as in this gauge perturbation theory has 
not been developed so far to the extent it has in other gauges: In the Schwarzschild case, most analyses of metric 
perturbations have been formulated so far within the Regge- Wheeler gauge |1C Jll|] (see, however, the recent mode 
decomposition of Schwarzschild's metric perturbations in the harmonic gauge ||].) In the Kerr case, so far the only 
practical approach for calculating the (mode-decomposed) metric perturbations is Chrzanowski's method |13|j , which 
is based on the radiation gauge. 

The main purpose of this paper is to provide a general prescription for calculating the gravitational self force in 
various gauges. To this end we shall first construct the general transformation law describing the behavior of the 
self force under a gauge transformation. Based on this transformation law, we re-express MSTQW's result ([[]) in an 
arbitrary gauge. We then re-formulate our method of mode sum regularization for a general gauge. 

The transformation rule describing the gauge transformation of the self force guarantees that the self force will be 
well-defined if (i) it was regular in the original gauge, and (ii) the gauge transformation is sufficiently regular (namely, 
the displacement vector £ M is sufficiently regular at the particle's location). A priori there is no guarantee that the 
transformation from the harmonic gauge to another desired gauge will satisfy this regularity criterion. One of the 
objectives of this paper is to explore the regularity of the self force in two commonly-used gauges: the Regge- Wheeler 
gauge and the radiation gauge. We find that the gauge transformations from the harmonic to these two gauges do not 
satisfy the required regularity criterion. As a consequence, our general transformation law does not yield a definite 
expression for the self force in these two gauges (the exception is the situation of a radial orbit in a Schwarzschild 
background, in which case the Regge- Wheeler self force is well defined). We note that this irregularity of the gauge 
transformation has been noticed independently by Mino [ |12| . 

This paper is arranged as follows. We start in Sec. || by exploring the way the gravitational self force transforms 



under a general gauge transformation. In Sec. Ill, which is somewhat out of the main course of our discussion, we 
consider the gauge transformation of linear gravitational forces in general. We find that this transformation law 
conforms with that of the gravitational self force. The general self-force transformation law is then used in Sec. [iy| to 
generalize MSTQW's expression for the regularized self force from the harmonic gauge to an arbitrary gauge. We also 
re-formulate our method of mode sum regularization for a general gauge. A few examples are provided in Sec. M, where 
we consider the transformation of the self force from the harmonic gauge to the Regge- Wheeler and to the radiation 
gauges. We find that in the Schwarzschild background the Regge- Wheeler self force is well-defined for a radial orbit, 
but is ill-defined for non-radial orbits. The situation with the radiation gauge is even worse: It is ill-defined even for a 
static test particle in flat space, and hence presumably also in all types of orbit in Schwarzschild or Kerr spacetimes. 



Finally in Sec. VI we summarize our main results and conclusions. We also discuss the indefinitcness of the self force 



in the Regge- Wheeler and radiation gauges, and suggest preliminary ways to overcome this difficulty. 
Throughout this paper we use metric signature ( — \- -\ — h) and geometrized units G = c = 1. 

II. GAUGE TRANSFORMATION OF THE SELF FORCE 

Our first goal in this section is to clarify the origin of the gauge dependence of the gravitational self force. Once 
this origin is well understood, the derivation of transformation law for the self force becomes rather straightforward. 

In discussing the origin of the gauge dependence, we find it useful to take the following point of view towards the 
gravitational self force kinematics: A point-like particle moves on a background metric go (e.g. the Schwarzschild 
geometry), and we wish to describe the particle's orbit. The particle, having a mass to, deforms the geometry, which 
is now described by the new metric, g = g + h, where h denotes the linearized metric perturbation produced by 
the particle. We also know that generally the particle will not follow a geodesic of go, due to its finite mass m. 
Since no external force is assumed to be present, one might attempt the simple point of view, according to which the 
particle moves on a "geodesic of the perturbed metric g" . This naive formulation, however, is unsatisfactory (if not 
totally meaningless), because the perturbed metric g is singular at the particle's location. We therefore must apply a 
different framework for analyzing the particle's motion: Assume that on the perturbed spacetime the particle follows 
a worldline x^(X), where A is an arbitrary monotonous parameter (we do not assume that A is a proper time in g, 
because the latter is not defined, due to the divergence of h). We now project the worldline x M (A) onto the background 
metric go on the basis of "same coordinate values" [we presume here that a choice of a coordinate system has been 
made in advance in each of the two spacetimes. Furthermore, we assume that the coordinates in the two spacetimes 
are "the same" if the small perturbation is ignored — which is equivalent to assuming that h is small, i.e. O(m)]. The 
projection defines a worldline x fJ, (X) on the background metric go, and we denote by r the proper time along this 
worldline (with respect to the metric go)- This construction now provides us with a natural definition of the self force: 
It is simply given by the acceleration associated with the worldline x M (r) in go, through Newton's second law: 

fd 2 x a „ , ,dx» dx"\ 

^ sm {-^ +Ta ^ww)- (2) 

In this expression, the connection T (just like the proper time r) is taken with respect to the background metric go- 
The origin of the gauge dependence of the self force is now obvious: Since g and go represent different geometries, 
in principle there is no unique way to project a point (or a worldline) from g to go- In the above formulation — as well 
as throughout this work — we adopt the rule of "same coordinate values" . Suppose now that an infinitesimal gauge 
transformation is carried out in the perturbed geometry g, associated with an infinitesimal displacement vector £**: 

x" -> x ,fi = x fl -^ (3) 

[this transformation changes h (and hence g), but of course the metric go of the background spacetime is unaffected]. 
The particle's worldline in the perturbed spacetime now takes a new coordinate value, ^(A) = x^(X) — £''. Projecting 
now the worldline on go, one obtains a new orbit x'^t'), where r' is the proper time (in go) of the new orbit 
x' M (A). It should be emphasized that the two projected worldlines, x^{t) and x'^{t'), represent two physically- distinct 
trajectories in go-Fl In particular, the self force will now take a new value, 

« = .»(^ + r;i„(x',——), (4) 

where T" (x 1 ) denotes the value of the connection in the new particle's location x la . 

We wish to calculate the quantity 5F^, which is the change in F^ l{ induced by the gauge transformation, to order 
to 2 [recalling that F s " lf itself is of order to 2 , and £ M is 0(m)]. To this end, we first transform the differentiation 
variable in Eq. (0) from r' to r: 



2 Recall, however, that in the perturbed spacetime g the two worldlines a: M (A) and x ,p, (X) are physically equivalent — they 
represent the same physical trajectory in two different gauges. This difference in the relation between :r M and x'^ in the two 
spacetimes simply reflects the non-uniqueness of the projection from g to go (which, in our "same coordinate value" formulation, 
is tied to the arbitrariness in choosing the gauge for h). 
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where 



dx' a /dr. Recalling that the term in squared brackets is already 0(m), we may omit the factor (dr/V') 2 — 



1 + 0(m), so at the required order we have 



f': 
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where (3 = m-^-^. Now, the force F s '" lf must be normal to the worldline (i.e., F' s " lf u' a — 0) by its definition in Eq. (0). 
We can therefore calculate it by projecting our last result on the direction normal to the worldline. Noting that the 
term j3u' a contributes nothing to this projection, we obtain 



F scu = m{6 x +u u x ) 



d 2 x' x 
dr 2 



+ T^(x') 
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dr dr 



Rewriting F^ u in the same form but with all primes omitted, and subtracting it from Fj^fe [evaluated at x'{x)\, we 
find at order m 2 



SF s a cU = m(5 a x + u a u x ) {q' x - q 



where 



,/A _ 



d 2 x' x 

~d7 2 ~ 



+ L ^ X) dr dr 



and q x is the same but with all primes omitted. [The term proportional to u' a u' x — u a u x does not contribute at the 
relevant order, because it is itself proportional to £ M , and q x and q' x are both 0(m).} All we now need is to calculate 
q' x — q x to leading order in ^ [expanding T x (x') about a; M to leading order in £**]. This is a standard calculation (it 
is often done when constructing the Jacobi equation for geodesic deviation), and one finds 



'A 



<l 



R\ a ,u^ a u l 



where an overdot denotes a covariant differentiation with respect to r and R x Ma „ is the Riemann tensor associated 
with the background metricn Now, the term u a u x in the above projection operator yields vanishing contribution 
when applied to the term including the Riemann tensor, due to the antisymmetry of the latter. Therefore, the final 
result is 



<50 = -m 



[g +u u )^\ + R ^xuu^t, u 



(6) 



(Since the calculation is carried out here at order m 2 only, in the last expression we may replace g by <7o-) 

The important message that arises from our discussion so far, is that the gravitational self force is a gauge-dependent 
notion. Specifying F^, u (t) by itself tells us almost nothing about the physical self force. In order for the information 
on the self force to have physical meaning, one must accompany it by the information on the gauge in which F s " lf 
was derived. Putting it in other words: The meaningful description of the gravitational self force must include both 
F S glf and the metric perturbation h a p. (Obviously, h a p contains the full information about the gauge.) This is closely 
related to a more general feature of general-relativistic kinematics (in the non-perturbative framework): Specifying 
the coordinate value of a worldline x^(t) tells one almost nothing about the physical nature of this trajectory, unless 
one is also given the metric g a p associated with the coordinates x^ . 

A remark should be made here concerning the regularity of the gravitational self force in various gauges. The 
construction by MSTQW yields a regular, well-defined, self force in the harmonic gauge. Therefore, in a given gauge 
G, the self force will be well defined if and only if SF^ lf is well defined. Obviously, if the gauge transformation from 
H to G is defined through a perfectly regular vector field £ A , the force in the G-gauge will be well defined. In most 
commonly used gauges, however, the vector field £ A associated with the transformation from the H-gauge to the 



J We use here the convention of Ref. ]14] for the Riemann tensor. Notice the different convention used by Mino et al. in 111. 



G-gauge may inherit some of the irregularity that h H itself possesses at the particle's location (to an extent that may 
depend on the gauge G and on the physical situation). In Sec. [V] this situation will be demonstrated for the RW 
gauge and for the radiation gauge. 

A priori it is not completely obvious what degree of regularity must be imposed on £ A in order for the self force to 
be regarded "regular" . Equation (ra) suggests a natural criterion for regularity: One should demand that £ A will be 
well defined (i.e. continuous) on the particle's worldline, and, furthermore, that along the worldline £ A will be a C 2 
function of r. Note, however, that there is some arbitrariness in choosing the regularity criterion. For example, one 
might impose a stronger regularity criterion, which requires £ A to be a C 2 function of x^ (such that the change in 
the connection due to the gauge transformation will be well defined); but we do not see much justification for such a 
strong demand. On the other hand, one may ease the above regularity criterion by extending the standard MSTQW 
regularization procedure and adding to it the element of averaging the self force (at a given moment) over all spatial 
directions. With this extended procedure of regularization, one may relax the demand for continuity of £ A at the 
worldline, replacing it by the weaker requirement that at the particle's location £ A will have a continuous limit along 
each spatial geodesic intersecting the worldline, and that this directional limit will be integrable over the solid angle. 
We further discuss this possibility at the end of the paper. 

For concreteness, throughout the rest of this paper we shall adopt the criterion which naturally follows from 
Eq. (ph — namely, that £ A be continuous on the particle's worldline. The second half of this criterion — the smooth 
dependence on r — will automatically follow, provided that the background metric (and hence also the particle's 
geodesic) is sufficiently smooth, which we assume here.[| 

III. GENERAL GRAVITATIONAL FORCES AND THEIR GAUGE TRANSFORMATION 

The above result (|6|) provides the full prescription for gauge-transforming the gravitational self force. It will be 
instructive, however, to address this issue of gauge transformation from yet another point of view, by introducing the 
notion of a (linearized) gravitational force and studying how this force transforms in a general gauge transformation. 

Consider again a spacetime described by a metric g = go + h, where go is a given background metric and h denotes 
a linearized metric perturbation. We do not assume in this section that h is a perturbation produced by a point 
particle; rather, h is assumed to be a prescribed weak gravitational perturbation (it may represent, for example, an 
incident gravitational wave). Suppose that a test particle with a mass m is moving freely in the perturbed spacetime. 
Obviously, this particle will move along a geodesic of g (we neglect the self force throughout this sectionPj). Namely, 
we shall have, in a given coordinate system x a , 

cPx* , a dx» dx v _ 

dr' 2 + "" dr' dr' ' { ' 

where x Q (r') denotes the particle's trajectory in the perturbed spacetime, r' is an affine parameter (with respect to g) 
along that trajectory, and T " are the connection coefficients associated with the metric g. However, we now wish to 
take the point of view according to which the particle traces a trajectory on the background metric go. This trajectory 
will deviate from a geodesic of the background metric go, and we shall interpret this deviation as representing an 
external "gravitational force" F" v , exerted on the particle by the perturbation h a p. This (fictitious) gravitational 
force is naturally defined as 

dx^ dx v 
d ,-'- » v ~dT~<fr- 



F« av = mx a =m(- j - r + 1^ — — ) , (8) 



where r is an affine parameter in the background metric go, an overdot denotes covariant differentiation (in g ) with 
respect to r, and r£o are the connection coefficients associated with the metric go- We wish to calculate F™ v to the 
first order in h (and to the leading order in to). 



4 We exclude here the situation in which the gauge condition defining the G-gauge explicitly depends on x M or r, and this 
explicit dependence artificially introduces non-smoothness to £ a (t) . In such spurious situations we must explicitly demand 
that £ A (r) be C 2 . 

5 Throughout this section we carry out the calculation to first order in the prescribed metric perturbation h, and to 
order in m [e.g., order m° in Eq. (M) below], so the self force is not included. 



A remark should be made here concerning the relation between the gravitational self force and the fictitious external 
gravitational force considered here. Obviously, the two notions are closely related, as both are defined through a 
mapping of a worldline from the physical spacetime g to a background metric go- Both forces are proportional to m 
and to the metric perturbation h (though in the self-force case one assumes that h is the metric perturbation produced 
by the particle itself). One may therefore be tempted to regard the self force as a special case of the more general, 
linearized gravitational force defined here. This is not quite the case, however. The gravitational force considered 
here is, after all, a fictitious force; that is, the particle actually follows a geodesic of the true physical metric g . This 
cannot be said about the orbit of a particle moving under the influence of its own gravitational self force: Since the 
self perturbation h is singular at the particle's location, the statement that the particle follows a geodesic of g = go + h 
is physically meaningless .£] For this reason, we must view the gravitational self force as a genuine, non-fictitious, force 
(though a delicate one, as expressed by its being gauge dependent). 

Proceeding with the calculation of F° v , wc first transform the differentiation variable in Eq. ([?]) from r' to r 
[mathematically this operation is the same one applied in the previous section, Eq. (0), though here it has a somewhat 
different meaning]. We find 

d 2 x a ^dx^dx" I ' dr 



dr 2 ' ^ dr dr + \dr ) dr' 2 dr °" ^ 

Denoting AT" = T" — T" and u a = dx a /dr, and substituting Eq. (0) in Eq. M) (keeping only terms linear in h), 
we obtain 

i^ Q rav = -mAr>"/ - Pu a . (10) 

We now get rid of the term f3u a by projecting FEL m on the subspace normal to u a , in the same way we treated SF^, lf 
above (recalling, again, that by definition F" av is normal to u a ). This yields 

*£.v = -m{5 a x + u a u x )AT^u^. (11) 

Expressing Ar in terms of h, we finally find 



F gr, v = ~~™(9 aX + u a u x ) {hxw + h Xv ;„ - V;a) U"""- ( 12 ) 



This expression (like the similar expressions below) is valid to linear order in the perturbation h, and on its right- 
hand side we may replace g aX by g$ A - It may also be useful to express Eq. (|l^) in terms of the trace-reversed metric 
perturbation h a p = h a p — ^g a ph (where h = g a ^h a ^). One easily obtains 

F g Q rav = mk a ^% r , s , (13) 

where k al3 ^ is a tensor given by 

k a ^ 5 = -g^u^u 1 - g^vCu 5 - ^u^u^u 6 + ^-u a g f ^u s + ^g aS g 01 . (14) 

Next we investigate how this gravitational force is modified by a general gauge transformation (ph. The metric 
perturbation h transforms according to 

h a p — > h a/3 = h a p + 5h a p , 

where 

5h a j3 = £ Q;/ 3 + £/3 ;Q . (15) 



6 One may take the point of view that the orbit of a particle under its gravitational self force is a geodesic in a spacetime 
with a metric go + /itaii, where /itaii denotes the tail part of the metric perturbation. This is, however, a fictitious geodesic, 
because the actual metric is go + h, not go + fttaii- (Recall also that in general h ta n fails to be a vacuum solution of the linearized 
Einstein equations.) 



From Eq. (|12|), the change in h will induce a corresponding change in the gravitational force F" v , given by 

5F£ av = -\m{g aX + u a u x ) {5h X)Jt , v + 6h x ^ - 5h^, x ) u»u v . (16) 

Do the self force F" lf and the linearized gravitational force F" av transform in the same manner? Substituting Eq. 
( |l5| ) for Shap in Eq. ( |lq ) and using the anti-commutation relation £, kXi , — £,p,;u\ = £,pR p p,\u, one obtains 

8F£„ = -m(g aX + u a u x ) (& i/u , + Z p R% Xv ) u»u» 

= -m [(g aX + u a u x ) 6 + K* M xvu"Z X u u ] . (17) 

Comparing this expression to Eq. (@) , we find that the two forces admit the same transformation law: 

SF^ V = SF s a clt . (18) 

This result is not surprising, because the two types of forces share a common kinematic feature: They are both 
constructed through a projection of a worldline from a physical metric g to a background metric go, and therefore 
they transform in the same manner. 

IV. REGULARIZING THE GRAVITATIONAL SELF FORCE IN VARIOUS GAUGES 

The method developed by MSTQW for regularizing the gravitational self force is formulated within the framework 
of the harmonic gauge. This means that in Eq. (fy) above, the two quantities on the right-hand side, i^ Q ( barc ) and 
F a ( mst ) are to be evaluated in the harmonic gauge — and the outcome is the self force in the harmonic gauge. We 
therefore rewrite this equation explicitly as 

p( H ) _ p( H ) _ p( H ) (to) 

^self — "'bare Mnst > \ 1J ) 

where the parenthetical index "H" denotes the harmonic gauge (for brevity we omit the tensorial index a here and 
in the equations below). 

Assume now that a gauge transformation is made, from the harmonic gauge to a new gauge which we denote 
schematically by "G" . According to the discussion in Sec. 0, the self force in the new gauge will be given by 



F^ = F™ + SF^r G) = 

self self ' self 



SF (H ^ G) + F (H) 



F ( H ) 
inst ' 



where SF^ el f* is the expression given in Eq. (p), with £ A being the displacement vector that transforms from the 
harmonic gauge to the new gauge G. To evaluate the term in squared brackets, we first recall that the "bare force" 
is related to the trace-reversed metric perturbation through 

Karc = mk a ^% T , S (20) 

[see the second equality in Eq. (28) of Ref. ||], which is expressed in terms of the metric perturbation itself as 

Karc = -\™{9 aX + u a u x ) {h X ^ + h x ^ - V;a) «"«"• (21) 

Noting further that Eqs. ( |l6| ) and (|lj) imply 

6F« n = -^m(g aX + u a u x ) {5h Xfi .„ + 6h x ^ - 5h^, x ) u"u v , (22) 

we then obtain (using hP^ + <Jfc,( H_ * G ) = h^) 

Fl"l + ^r G) =FiZ (23) 

This result has a simple interpretation in terms of the notion of "gravitational force" discussed in the previous section: 
(i) As was established there, the self force and the gravitational force transform exactly in the same manner, and (ii) 



the "bare force" is nothing but the gravitational force associated with the full metric perturbation h (produced by 
the particle). It then follows that the self force and the bare force transform in the same manner .Q 
We conclude that in an arbitrary gauge G the regularized gravitational self force is simply given by 

p( G ) _ p( G ) _ p>( H ) (OA) 

Namely, in an arbitrary gauge G, the singular piece to be subtracted from the bare force is always the instantaneous 
piece expressed in the harmonic gauge, and not in the gauge G, as one might naively expect. 

Our last result is of special importance: The analysis by MSTQW tells us how to calculate the physical self force 
associated with the metric perturbation in the harmonic gauge. In particular, it tells us how to construct the "correct" 
instantaneous part of the bare force in this gauge. Our above discussion implies that even when calculating the self 
force in a different gauge, the "correct" instantaneous part must still be calculated in the harmonic gauge. [The 
explicit construction of the instantaneous part from the harmonic gauge Green's function is described in Eq. (29) of 
Ref. ||]. This harmonic-gauge-related instantaneous part is the one which captures the "correct" divergent piece to 
be removed from the bare force in whatever gauge. Intuitively, this special significance of the harmonic gauge may 
be attributed to its inherently isotropic nature: The "correct" divergent piece that should be removed from the bare 
force must be spatially isotropic (see, e.g., the analysis by Quinn and Wald &), and it is the harmonic gauge which 
admits this isotropic structure; other gauge conditions may introduce an artificial distortion to the singular piece. 

Mode-sum regularization in various gauges 

In Ref. p| we introduced a practical calculation scheme for the gravitational self force, based on the regularization 
procedure by MSTQW, which employs a multipole mode decomposition. This method of "mode sum regularization" 
has been formulated in Ref. fq] only within the harmonic gauge. Let us now examine how the above discussion, 
concerning the construction of the regularized self force in various gauges, applies in the framework of the mode-sum 
scheme. 

Within the mode sum scheme, the harmonic-gauge regularized gravitational self force is given by H 



F a(H) _ V- 
r sclf — / , 



Z ( F bir ( c H) - A " L -B a - C a /L) - D a , (25) 

=0 



al(B) 



where the summation is over multipole modes I, and L = 1 + 1/2. In this expression, F hs ^. c is the contribution to the 
self force from h} , the Z-mode of the metric perturbation in the harmonic gauge. This contribution is given by 

<£ H) = mk^h l ^ 5 , (26) 

where ft.i is the trace-reversed /ii and k af3l& is the tensor given in Eq. (|l4|). The vectorial quantities A a , B a , 
C a , and D a appearing in Eq. ( |25| ) are /-independent. These quantities, which we call "regularization parameters", 
are constructed from the Z-modes of the instantaneous part F-lJ , in a manner described in Ref. M . 



The prescription provided by Eq. (25) yields the "harmonic gauge" self force. It is now possible, however, to 



re-formulate this prescription in any other gauge "G" , using 

F S T=F S T + 6F^\ (27) 

Rewriting Eq. ( |22] ) as 5F£, lf — mk a P jS 5hp 1] s (where Sh denotes the trace-reversed Sh), and decomposing Sh into 
/-modes, we obtain 

oo 



7 It should be emphasized that the physical notion of "gravitational force" introduced in the previous section is not necessary 
for the derivation of Eq. (|23|). Thus, starting from Eq. (pi), one can derive Eq. (G2J) directly as a mathematical identity 
[following the same mathematical steps used above for constructing Eq. (jlTf) from Eq. fllq)], without any reference to the 
notion of "gravitational forces". 



Substituting this and Eq. ( |25| ) into Eq. fl27|), we find 

oc 

F s Q j f G) = J^ [(<r ( c H) + mk a ^ s 5h l ^ G) ) - A a L -B a - C a /L~\ - D a . (28) 

1=0 

Using now Eq. (|26|), we can re-express the term in parentheses as 

mfe^^S + mk^ s Sh l ^ G) = mk°*"h l ™ = F*£\ 

where Ff* J^ denotes /-mode contribution to the "G-gauge bare force" , namely, the contribution to the bare force 
from the mode I of the (bare) metric perturbation in the G-gauge, through Eq. (£0|). We thus obtain the simple 
expression for the self force in an arbitrary gauge "G" , 



self — 

1=0 



oo 

J2 ( F bl ( c G) ^ AaL -B a - C a /L\ - D a . (29) 



We conclude that the regularization parameters A a , B a , C a , and D a are independent of the gauge. This result 
has a simple intuitive explanation: These parameters are determined by the mode decomposition of the instantaneous 
piece of the metric perturbation, which — based on our above discussion — is always to be expressed in the harmonic 
gauge, regardless of the gauge chosen for calculating the self force. Thus, the regularization parameters A a , B a , C a , 
and D a are, in effect, gauge-independent. 

It should be commented that the above discussion is valid as long as 5F ' "L ~* (and hence the self force in the 
gauge "G") admits a well defined finite value. As we demonstrate in the next section, in certain gauges the self force 
turns out to be irregular or ill-defined. In such cases, the irregularity may enter Eq. ( |29|) through the bare modes 
■^bare anc ^ render the sum over I non-convergent. 

V. EXAMPLES 

In this section we study the transformation of the self force from the harmonic gauge to other, commonly used gauges, 
in a few simple cases. In principle, this transformation is done by first solving Eq. (FFq) for the gauge displacement 
vector £ M , and then constructing the force difference <5F s " lf by using Eq. (J6). We shall primarily be concerned here 
about the regularity of the self force in the new gauge. As discussed in Sec. H], we shall regard the G-gauge self force 
as regular if the vector field £ p is continuous at the particle's location. If it is indeed continuous, then the self force 
in the new gauge is given in Eq. (o) (the demand for a C 2 dependence on r is automatically satisfied, as discussed in 
Sec. [n]). We begin by considering the transformation to the Regge- Wheeler (RW) gauge, for radial trajectories in the 
Schwarzschild spacetime. Then we examine the transformation to the RW gauge for a uniform circular orbit. Finally, 
we examine the transformation to the (outgoing) radiation gauge, in a simple flat-space example. 

A. Regge— Wheeler gauge: radial trajectories 

We consider a particle of mass m moving along a strictly radial free-fall orbit on the background of a Schwarzschild 
black hole with mass M 3> m. (Of course, the motion of the particle will remain radial even under the effect of 
self- force, by virtue of the symmetry of the problem.) In what follows we use Schwarzschild coordinates t, r, 6, ip and 
assume, without loss of generality, that the radial trajectory lies along the polar axis, i.e., at 9 = 0. 

Let fv a J and h^* denote the metric perturbation produced by the above particle in the harmonic and RW gauges, 

respectively. The displacement vector field ^ which transforms h a J to h a * satisfies the gauge transformation 
equation 

C w) = C )+ ^+^- ( 3 °) 

The symmetry of the physical setup motivates one to consider only axially-symmetric even-parity metric perturbation 



modes. Accordingly, we shall look for solutions to Eqs. (30) which are ^-independent and also have £ y = 



For even-parity perturbation modes, the RW gauge conditions take the simple algebraic form |1( 



, (RW) _ , (RW) _ , (RW) _ n 



(31) 



where /i ang = (hgg — sin 2 $h vv ) /2. Imposing these conditions, the gauge transformation equation ( |30| ) yields three 
coupled differential equations for the three components £ t , £ r , and £#: 



6,0 + 6 



,(H) 

he ' 



(32a) 



&-,fl+6,r-(2/r)& = -ft 



(Hi 



(32b) 



sin (sin 1 0^e 



"ang ■ 



(32c) 



Eq. (32c) can be immediately integrated with respect to 9 (with fixed t,r), yielding 



& 



sin 6* 



sin- 1 d'h^Ue' 

ang 



Mr,t) 



(33) 



where ?/>i is an arbitrary function. (As we shall discuss belo w, h an l va nishes sufficiently fast as 9' — ► 0, such that the 
integral is well-defined at the lower limit.) Then, Eqs. (32a) and (B2q) are immediately solvable, yielding 



6 = - 



i (H) 

he 



+ &,t dO 1 + M^ t), &■ = - / Kg' + & r - (2/r)£,j d0' + M^ t), 



,W 



(34) 



where "02 and ^3 are two other arbitrary functions.^] 

Now, in order to explore the behavior of the quantity SF" lt corresponding to the gauge transformation H^ RW, 
by means of Eq. (g), one has to characterize the behavior of the vector field £ M at the location of the particle. This 
requires one to first explore the behavior of the various H-gauge metric functions appearing in Eqs. p3) and (pi) 
at the particle's location. This task is most easily accomplished by considering the Hadamard form of the metric 
perturbation in the neighborhood of the particle. For the trace-reversed metric perturbation in the harmonic gauge, 
this form was given by Mino et al. [sec Eq. (2.27) of Ref. ||; Alternatively, see Eq. (45) of Ref. (|]]:| 



:(H) 

'a/3 



4me 



*Uf3 



0(e°), 



(35) 



where e is the spatial geodesic distance to the particle's worldline (i.e., the proper length of the geodesic normal to 
the worldline which connects the latter to the evaluation point), and the terms included in O(e ) are assured to be at 
least C 1 functions of the coordinates at e = 0. The metric perturbation itself is then given by 



,(H) 
l af3 



4me * {u a U(3 + g a p/2) + O(e ). 



(36) 



Since the worldline is radial (namely ug 
and /ling appearing in Eqs. ( p3| ) and * 

i\(~\rr\ oil romilcir li ^i of loacf f^*-\ r\n 4- 



'he> 



u v = 0), it now follows that the metric perturbation components 

hrfj' , and ft^Sg appearing in Eqs. ( |33| ) and ( |34| ) all have vanishing contributions from the singular 0(e _1 ) term, and are 
therefore all regular (i.e., at least C^) on the worldline. Consequently, one can easily construct solutions for f , which 
have regular, finite values at the particle's location: Starting from Eq. (33), we first observe (e.g., by transforming 
to cartesian- like coordinates at the polar axis, and demanding axial symmetry as well as C 1 asymptotic behavior at 
9 = 0) that /ling falls off at 9 — > faster than 9. As a consequence, the integral in Eq. (p3|), too, falls off faster than 
9. Thus, £e is regular at 9 — ► 0, and it vanishes there like ex 9. (With the choice ipi = 0, £g would vanish even faster 
than6> 2 ). 



8 The arbitrary functions ipi represent a true freedom in the construction of the RW-gauge metric perturbations. This may 
be attributed to the freedom of specifying the monopole and dipole modes of the metric perturbation — see the discussion in 
Ref. |§. 

9 To obtain Eq. (B5I) from Eq. (2.27) of Ref. |l|, recall that at the location of the particle we have g 12 a = 8£ and k — 1 (using 
the notation of M). 
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Consider next the two integrals in Eq. (|3J) . From the above discussion it immediately follows that the two derivatives 
£$ t t and £,e,r vanish like oc 9 (at least) — like £g itself. Since hf e and hf e are regular (C 1 ) too, we find that the two 
integrands in Eq. ( |34| ) are bounded at 6 = 0. (In fact, by transforming to cartesian-like coordinates near 9 = one 
can easily verify that h^ e and h^g — and hence the two integrands — vanish at 6 — ► 0.) Consequently, the two integrals 
vanish at 9 — > 0. We find that along the particle's worldline all components of £ M are regular, and satisfy 

£,e = 0, & = ip2(r,t), £ r = ipa(r,f), 

where fair, t) and ip3(r, t) are freely-specifiable functions. (In fact, this holds not only at the particle's worldline, but 
everywhere along the polar axis.) Furthermore, choosing i\>i = "03 — 0; we obtain a solution for £** which is not only 
regular but is also vanishing along the particle's worldline: £ M (r) = 0. 

Since the above-constructed vector £ M is continuous at the particle's location, we obtain — through Eq. (|g) — a 
regular finite value for the desired quantity 8F^. Thus, for strictly radial trajectories in Schwarzschild spacetime, 
the gravitational self force is regular in the RW gauge. Moreover, this RW-gauge self force can be made equal to 
the harmonic-gauge self force, by exploiting the remaining freedom in the RW gauge (manifested here by the three 
arbitrary functions ^1-3). 

B. Regge— Wheeler gauge: circular orbits 

Let us now consider a particle which (in the lack of self force) moves on a circular geodesic at r = r^ > 6A/ around 
a Schwarzschild black hole. Without loss of generality, we shall assume an equatorial orbit (i.e., 9 = ir/2 and u e = 0) 
and will consider the self force at a point P located on the particle's orbit at t = <p = 0. In this physical scenario, 
the metric perturbation contains both even and odd parity modes. The RW gauge condition |1C| ] then becomes a 
bit more complicate than the one specified in Eq. ( pl| ) for a purely even perturbation (in general, the two algebraic 

conditions h t8 = h r9 = are no longer valid, and are to be replaced by conditions involving derivatives of the 
metric perturbation). However, the two gauge conditions involving the angular components of the metric perturbation 
maintain a simple algebraic form, namely 

C° = °. hi ™ } - °- ( 37 ) 

For our purpose, it will be sufficient to consider only these two conditions. When imposed on the gauge transformation 
equation ([30|), these conditions lead to a set of two coupled equations for £g and £ y : 

sin0(sin- 1 6»ee),e -sin - ^^ = -ftW. (38 a ) 

Ze, v + sm 2 9(sm~ 2 9Z v ),e = -h < ff, (38b) 

The source terms for these equations are evaluated, again, with the help of Eq. (pq): We find that h e is regular 



at the particle's location, but the source for Eq. (38a) diverges there as — h an g — ar 2 e , where a = 2mr ( 



^ ' r- ° t-i rv^ r J- TrrVt rvf ir\ /-i Oivn/r< II 



2m(ro/M — 3) _1 Q and, as before, e denotes the spatial geodesic distance to the particle's worldline. In what follows 
we analyze the behavior of £g and £ v at the immediate neighborhood of P, to leading order in e. 

We first note that no derivatives with respect to r and t appear in Eq. ( |38a| ) (though the source term depends on 
r and t through e). Therefore, this equation can be solved for each r,t separately. For our purpose — demonstrating 
the discontinuity of the solution at P — it will be sufficient to consider the solution at the two-dimensional plain 
r = fo, t = 0, which is simpler to analyze. 

To bring Eqs. ( p8|) to a convenient form, we introduce the local cartesian-like coordinates y = ro sin^sintp, z = 
ro cos 9 in the neighborhood of P. Note that z = y = at P, and that (for r = ro,t = 0) at the leading order we 
have e = [(1 — v 2 )~ 1 y 2 + z 2 ] 1 / 2 . Here v denotes the particle's velocity in the Lorentz frame of a static local observer, 
v = (—g V{p /gtt) l l 2 {d(p/dt). One can easily obtain the explicit value of v. |l7j 

» = (-9 vv /9 U ) 1/2 (u v /ut) = (r /M - 2)- 1 ' 2 < 1. 

Transforming in Eqs. (|38| ) from (9, ip) to (z, y) we obtain two coupled equations for £ z and £ y , reading 

(1 - z 2 /rl)^ z - (1 - y 2 /r 2 )Z y , y = a/e + ■ ■ ■ , (39a) 
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sz,i/ + Sy,. 



2yz{rl - z 2 ) l £ VfV 







(39b) 



where the dots (• • •) represent corrections to the source term which are at least C 1 at P. As we are interested only 
in the leading-order behavior of £ M at P [where (z/ro) 2 , (y/Vo) 2 , (zy/r 2 ) all vanish], we shall proceed by restricting 
attention to the leading-order form of Eqs. (39): 



Kz.z t,y 



a/e, 



(40a) 



Sz,y ~r ^y,z U. 

Equation (40b) allows us to express the vector £ M in terms of a scalar potential <&, as^j 



With Eq. (40a), this potential is then found to satisfy Poisson's equation 



$ 



$ 



a/e. 



(40b) 

(41) 
(42) 



It is convenient to introduce polar coordinates in the zy-plain, which we define through z — psin0, y — pcosi 
Transforming in Eq. (142) from z, y to p, (j), and substituting e — [(1 — v 2 )~ 1 y 2 + z 2 ] 1 / 2 , wc obtain 



P' 1 Q>*,p),p 



p- 2 ^. H = -(l-v 2 )^ 2 (l-v 2 sm 2 



-1/2 



(43) 



Next, we wish to expand <&(p, 4>) into angular Fourier modes e m< ^. Before doing this, however, there is a subtlety 
that must be discussed. The displacement vector £ M must be a single-valued (SV) function of </>. This means that 
both $ ? and <£ iP must be SV too. However, in principle the generating potential $ needs not be a SV function of 
(j). Therefore, in the complete mode decomposition of $ one may also include certain functions of 4> which are not 
necessarily SV. However, since the </>-derivative of each such multi-valued function must be SV, this function must 
be linear in (j) (such that the Fourier expansion of 3?^ will only include SV Fourier modes). Furthermore, since the 
p-derivative must be SV too, this "linear mode" must be independent of p. The full decomposition thus takes the 
form 



$(p,0)=«/>+ J2 e m ^n{p), 



(44) 



where c is an arbitrary constant. Substituting this form in Eq. (43) (recalling that c(j) satisfies the homogeneous part 
of this equation), one obtains an ordinary equation for each n-mode, reading 



P Up^r. 



- n 2 p 2 $„ = - /„ , 



(45) 



where the coefficients /„ are given by the (elliptic) integrals 

,-2tt 



In 



VT 



2n 



-iri(p 



v/T 



v 2 sin i 



(46) 



It can be easily verified that /„ vanishes for all odd n. For even n, however, /„ is generally non- vanishing. In 
particular, for n = the integrand in Eq. (Eq) is bounded from below by unity, hence fo > vl — v 2 > 0. 
The general exact solution to Eq. (Esl) is easily constructed: 



$ = 



b p + a +/3 ln/9, 



for n = 0, 



b n p + a nP W+l3 n p-W, foin^O, 
where a n and (3 n are arbitrary constants, and 



(47) 



10 Defining E = (E y ,E z ) = (— £ y ,£ z ), Eq. (40b) reads V x E = 0, which allows one to define E = V3>. (Note, however, that 
since there is a singularity at y = z — 0, $ needs not be single- valued - see the discussion below). 
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a/n/(l — n ) for even n, 
for odd n. 



We may now construct the modes of £ M by applying Eq. ([41J) to each of the single modes. We then wish to figure out 
what is the solution with the most regular behavior at the limit p — > 0, which concerns as here. Clearly, any choice of 
[3 n 7^ will lead to a divergent $ n and hence to a divergent vector £ M . [Note that the norm of (£ w ,£z) is the same as 
that of V<& , and is hence bounded below by |3?, p |; and the contribution to the latter from a nonvanishing (3 n would 
diverge like ex p - '™' -1 .] Similarly, a nonvanishing c would yield a potential $ whose (normalized) derivative in the 
tangential direction, p _1 $^, diverges like cp _1 .Fj The most regular solution is thus one with /?„ = for all n, as well 
as c = 0. This solution takes the form 

fc(p,0) =a + pH(<f>)+0(p 2 ), 

where 



H(<t>) = (ale 2 * + a_ ie -**) + ]T 6„ 



ira<£ 



Returning from <I> to £ M , we find e.g. for the cartesian-like component S, y (ignoring higher-order contributions in p): 

€v = ~®,v = -P,v H - P<t>,yH,4, . 

Substituting p y — costfi and (f). y = —p^ 1 sin</>, we find 

£„ = -Hcos<p + H^sm<j) = £ y (<l>) . 

Clearly, in order for £ y to be continuous at p — ► (where (/> is indefinite), it must be independent of <p. However, 
i y ^ = {H^tf, + H) sin^, and 



3C 



E M + H= J2 0- - ™ 2 )KJ n * =a J2 ^ e 



intp 



n— — oo n— — oo 



This function of <j> does not vanish (identically) unless all coefficients /„ vanish; however, as was shown above, /q ^ 0. 
We find that £, y (p — > 0) does depend on <\> (the same can be shown for £ z ). This means that the vector £ 



/< 



discontinuous at P.rj 

As the gauge displacement vector £ M does not admit a definite value at the particle's location, Eq. (|6|) cannot be 
used, as it stands, for constructing the self force in the RW gauge. Following the discussion at the end of Sec. |J, 
we arrive at the conclusion that in the case of circular motion, the "RW self force" is ill defined (unless one further 
extends the regularization procedure — e.g., by introducing an average over solid angle; see the discussion in Sec. O). 

We conclude this discussion with two remarks: First, though the discontinuity of ^ was explicitly demonstrated 
here for circular orbits, this conclusion should also apply to generic non-radial, non-circular, orbits (for radial orbits, 
however, it was demonstrated above that £ M is continuous). 

Second, the above construction shows that for a suitable choice of the free parameters (namely c — (3 n = 0) the 
component S, y is bounded at P. The same holds for £ z . This implies that £g and £ v are bounded (though discontinuous) 
at the particle's location. It still remains to be checked, however, whether £ t and £ r are bounded or not. 



11 Divergent contributions from different n-modes cannot cancel each other, because they have different dependence on 0, as 
well as different rates of divergence (p~' n ' _1 ). Also, a divergence coming from the linear mode cannot cancel a divergent n = 
mode, even though in both modes |\7<3>| ex p _1 , because the direction of V$ is tangential for the linear mode and "radial" for 
the n — mode. 

12 The indefiniteness of the RW self force could be intuitively understood, by realizing that the RW gauge condition "distracts" , 
to some amount, the presumed isotropic structure of the divergent local piece of the metric perturbation, by artificially signifying 
the 6 direction. (This isotropic structure is best accounted for within the harmonic gauge.) 
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C. Radiation gauge 

Finally, we examine the transformation of the self force to the so-called "radiation gauge". (We recall that, so far, 
the mode decomposition of metric perturbations in Kerr spacetime has been formulated primarily within the radiation 
gauge |13j.) We shall consider here the simplest possible case: a static particle in flat spacetime. As we shall shortly 
see, even in this trivial case, the gauge transformation from the harmonic to the radiation gauge is pathological, and 
the metric perturbation (and hence the self force) is ill defined. 

We shall specifically consider the outgoing radiation gauge (similar results are obtained when considering the ingoing 
radiation gauge). We use standard flat-space spherical coordinates t,r,6,ip, and assume that the static particle is 
located off the origin of the spherical coordinates, i.e. at some r = tq > 0. Also, without loss of generality, we 
locate the particle at the polar axis, 9 = 0. The outgoing null vector field takes the form l a = [1, 1, 0, 0]. The metric 
perturbation in the radiation gauge, h a/3 , is defined by the requirement! 1 3 | 



Consider now the ^-component of Eq. m 



*3>v=o. 






which reads 






e } + >& 


= 0. 




= h a/3 + &*;/3 + 


£/3;a, 


this becomes 


; r>t + 2£ M = -h 


(H) 
tt 





(48) 



With the gauge transformation equation h a J 

&,r + £r,t + 26,*=-^ ->»£'■ (49) 

Motivated by the staticity of the problem, we shall only consider ^-independent solutions, so t- r t and £t,t may be 
dropped. Also, in the harmonic gauge we have 

4 H } = 2m/e, h™ = 0, 

where e denotes the spatial distance to the particle's location (this may be easily obtained by transforming the well 
known Coulomb- like Cartesian solution to spherical coordinates). Equation ( |49| ) now becomes 

&,r = -2m/e . (50) 

At this point we introduce standard Cartesian coordinates t, x, y, z , such that the particle is located at the origin (a; = 
y = z = 0), and the z direction coincides with the radial direction at the particle's location (namely, x = r sin#cos<p, 
y = rsmO simp, and z = rcosQ — tq). At the leading order in e, we may replace d r by the Cartesian derivative 
operator d z . Equation ( pOJ ) then becomes 

6, 2 = -2m(z 2 + p 2 )- 1 / 2 , (51) 

where p 2 = x 2 + y 2 . Eq. ( pT[ ) can now be easily integrated with respect to z (with x,y held fixed). We obtain 

6 = -2m log (- + ^l + z 2 /pA + R(x, y) , (52) 

where R(x,y) is an arbitrary function. This is the most general (t- independent) solution for ^ t - 
Consider next the asymptotic form of £( as we go to the limit x, y — ► with fixed z^0. One finds 

t (n^n\ c* / +^rn\n(p/2z) + R(x,y), z > 0, , . 

^ P U) -\-2mln(p/2\z\) + R(x,y), z<0. l ° dj 



13 In the case of a pure vacuum perturbation over a Kerr background, the additional condition h, = g a h a J = can be 
imposed in a consistent manner, as done by Chrzanowski in (l3(. Here we consider the perturbation in a region surrounding 
a point source, and it is unclear to us whether the additional condition hS ' — will be consistent with the gauge condition 
((tq). We shall therefore not make any use of this extra condition here. 
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By a suitable choice of the function R(x,y) one may, at best, eliminate the divergence along one of the rays z < 
or z > (by choosing R ~ ±2m\np, respectively), but not along both rays simultaneously. We thus arrive at the 
conclusion that £< unavoidably diverges logarithmically (at least) on approaching the axis p = 0, along either the 
z < ray or the z > ray (or both). Constructing now the tx and <y components of the radiation-gauge metric 
perturbation, we find h\ x = £t,x °c ^/p 2 ! an d a similar expression for /iL , as p — > (at cither z < or z > 0). 
Namely, the metric perturbation inevitably diverges at least along half the axis p = 0. 

It thus turns out that in the radiation gauge, the perturbation associated with a pointlike particle is represented 
by a string- like one-dimensional singularity. In particular, the radiation-gauge metric perturbation cannot be well 
defined in a complete neighborhood of the particle. (Compare with the harmonic or RW gauges, where the singularity 
is confined to the particle's location and the metric perturbation is well defined everywhere in the particle's neighbor- 
hood.) This pathological behavior — manifested already in the elementary case of a static particle in flat space — serves 
to demonstrate the pathological nature of the radiation gauge in the presence of point sources. As the radiation 
gauge seems inappropriate for representing the metric perturbation in the particle's neighborhood, it becomes rather 
meaningless to consider the self force acting on the particle in that gauge. 

Finally we note that although the indefiniteness of the radiation-reaction self force was demonstrated here only for a 
static particle in flat space, the same indefiniteness should also occur generically for all types of orbits in Schwarzschild 
or Kerr spacetimes. 

VI. SUMMARY AND DISCUSSION 



The main results of this manuscript are contained in Eqs. (|6|), (|24|), and (|29|). Eq. (p|) describes the gauge trans- 
formation of the gravitational self force, given the gauge displacement vector £ M . Eq. (p4J ) describes, in a schematic 
manner, the extension of the MSTQW formulation for the gravitational self force to an arbitrary gauge 'G': It implies 
that the "correct" singular piece to be removed from the bare force in the G-gauge [the one derived directly from 
the G-gauge metric perturbation through Eq. (Effl] is always to be calculated in the harmonic gauge, as described 
in the original analysis by MSTQW. By applying these results to our mode-sum regularization method (which was 
previously formulated only within the harmonic gauge la]) we finally obtained Eq. (p9|), which describes a practical 
mode-sum prescription for construction of the gravitational self force in any gauge 'G' (provided that the self force 
has a regular, finite value in that gauge). We stress again that, since the gravitational self force is a gauge-dependent 
notion, expressions like Eq. ( f24[ ) or Eq. (g9|) for the self force will be meaningful only when accompanied by the full 
information about the gauge to which they correspond. (Alternatively, one can specify the metric perturbation hS ' 
itself, which of course contains the full information about the gauge.) 

The implementation of Eq. (g9J) for calculating the G-gauge self force involves two distinct parts: (i) calculation of 
the bare modes of the force in the G-gauge [through Eq. (J20|)] ; and (ii) derivation of the four vectorial regularization 
parameters A a , B a , C a , and D a . Our discussion concerning the gauge transformation of the self force led us to 
conclude that the values of these regularization parameters do not depend on the gauge in which one calculates the 
self force: These parameters are always to be calculated in the harmonic gauge (using the analytic technique described 
in Ref. [pi). This "gauge invariance" property of the regularization parameters is demonstrated by the recent analysis 
by Lousto pjlq, who calculated (numerically) the values of A a , B a , and C a in the RW gauge, for a radial orbit on a 
Schwarzschild background. These numerical values appear to be in perfect agreement with the harmonic-gauge values 
derived analytically in Ref. || (in the case studied so far, of the self force at a turning point of a radial geodesic) . 
Also, the (zero) value obtained for the parameter D a in the harmonic gauge S agrees with Lousto's result for D a in 
the RW-gauge (which was based on a proposed zeta-function regularization procedure M ) . 



The prescription (29), as well as Eq. (|24|), is only applicable when the self force admits a definite finite value in 
the G-gauge. Whether or not this is the case for a given gauge "G", can be decided with the help of Eq. (0): The 
analysis by MSTQW implies that the self force will always have a regular finite value in the harmonic gauge (and it 
also tells us how to derive this value). Therefore, the G-gauge self force would be well defined, in our approach, only 
if the transformation from the harmonic gauge to the G-gauge would yield — through Eq. (|g) — a regular finite value 
for the force difference <5F s " lf . It is only in this case that we are able to use Eq. (G3) for calculating the G-gauge self 
force. Otherwise (namely, if SF^ U diverges or is indefinite), Eq. ( p9|) appears to be useless. 

As an example, in Sec. Mwe explored the transformation from the harmonic gauge to the Regge- Wheeler gauge. We 
found that the RW self force is well defined as long as strictly radial trajectories are considered. For such trajectories, 
Eq. ( p9| ) then provides a useful prescription for computing the RW self force. However, this seems not to be the 
case for more general orbits, as we demonstrated by considering a circular orbit: Here, the transformation from the 
harmonic gauge yielded an indefinite value for the RW self force. The situation is even worse in the radiation gauge, 
where SF^ { ~* is found to be not only discontinuous but also unbounded, and presumably for all types of orbits. 
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How could one interpret a situation where SF^ 1{ diverges (or is indefinite)? In some occasions, such a result may 
be attributed to a severe pathology of the gauge. This seems to be the case in the radiation gauge, as implied by the 
fact that in this gauge the metric perturbation diverges not only at the particle's location, but also along an (ingoing 
or outgoing) radial ray emerging from the particle (see Sec. |V|). However, the situation seems to be different in the 
RW gauge, in which the metric perturbation is well-defined in the neighborhood of the particle (though of course not 

at the particle itself), like in the harmonic gauge. In this case we have seen that, for non-radial orbits, SF"\ f ~* 

(and hence also F" lf itself) is ill defined. This originates from the fact that certain components of £ M — e.g. £ e or 
£, v — admit a direction-dependent limit (as demonstrated by the dependence of e.g. £ y on (j>; cf. Sec. |y|). 

This situation — a direction-dependent expression for the self force in certain gauges — motivates one to consider 
a simple generalization of the standard MSTQW regularization procedure, by averaging over all spatial directions. 
Namely, one can evaluate the limit of the right-hand side of Eq. ( p4| ) (or, similarly, the limit of the displacement 
vector £ M ) along fixed spatial null geodesies emanating from the particle, and then average over the solid angle (in 
the particle's rest frame). This would clearly be a generalization of the MSTQW procedure, because whenever the 
coincidence limit is well defined, the average over solid angle will be well-defined too, and will yield the same result. 
One still needs to investigate how this averaging over directions is to be implemented within the context of the 
mode-sum regularization. 

The above generalized regularization procedure will yield a definite self force in a wide class of gauges (though 
not in all gauges; Obviously one can construct a displacement vector £ M which does not even have a directional 
limit, in which case the generalized regularization procedure will fail to yield a definite self force). The analysis in 
Sec. |v| suggests that for circular orbits the displacement vector £ M from the harmonic to the RW gauge may have a 
well-defined directional limit, and hence the RW self force may be well defined within this generalized prescription. 
Recall, however, that the above analysis does not completely guarantee this regularity of the (generalized) RW self 
force, because so far we have only analyzed the tangential components £ y and £ 2 (which yield £ e and £ v ), but not 
£,t and £ r . Also, our analysis was restricted to the surface r = ro,t = 0, i.e. to directional limits through tangential 
directions. 

There seems to be another procedure that would allow one to use the metric perturbations in e.g. the RW or 
radiation gauges for useful self- force calculations (without resorting to the above generalized regularization procedure) . 
We shall now briefly outline here a preliminary version of this procedure. (We note that a similar approach has been 
proposed by Mino Q .) Suppose that the metric perturbation h^ is known (e.g. in the form of mode decomposition), 
where "G" refers to either the RW or radiation gauges. If we knew how to convert h^ G ' to the harmonic gauge, it 
would be straightforward to construct the self force from it, through Eqs. (119) or (pq). However, performing the 
transformation G— >H requires one to solve a system of partial differential equations for £ M , and unfortunately we do 
not know the exact solution of this system. Nevertheless, it appears possible to construct an approximate, leading- 
order, solution of this system, for both the RW and the radiation gauges. This was demonstrated in Sec. M (for both 
gauges) in a few simple cases, and it appears likely that the leading-order solution can be generalized to a generic orbit. 
Let us denote this leading-order solution by £^ . In principle one can then use £ M to transform the metric perturbations 
from the original gauge G to an "approximate harmonic" gauge, which we denote H. Presumably, in the gauge H the 
self force will be well defined, since the metric perturbations in the harmonic and H gauges share the same leading- 
order asymptotic behavior. After decomposing £ M into ^-modes, one can use the mode-sum regularization method 
to calculate the H-gauge self force: Applying Eq. ( p3[ ) for each of the single /-modes, with "H" and "G" replaced, 
correspondingly, by "G" and "H" , we first get 

pai(H) _ pai(G) Q /(G^H) 

-'bare — ^barc ' ox sclf ' 

where SF" el \ ~* is to be obtained from Eq. <|6|) by replacing £ A by the Z-mode of £ M . Then, writing Eq. ( f29[ ) for the 
H-gauge (i.e., with all "G" replaced by "H") and substituting the above expression for F^i , one obtains 



r solf — / , 



[(Xll^ + ^self^) " A * L ~ B ° ~ C ° '/ ' L ] ~ D" > ( 54 ) 



which provides a prescription for calculating the H-gauge self force through the modes of the bare force in the G 
gauge. We hope to further develop and implement this method elsewhere. 
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